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ABSTRACT 

In the present paper we construct maps of polarized synchrotron radio emission of a whole galaxy, 
based on local models of the cosmic ray (CR) driven dynamo. We perform numerical simulations of the 
dynamo in local Cartesian domains, with shear-periodic boundary conditions, placed at the different 
galactocentric radii. Those local solutions are concatenated together to construct the synchrotron 
images of the whole galaxy. The main aim of the paper is to compare the model results with the 
observed radio continuum emission from nearly edge-on spiral galaxy. 

On the basis of the modeled evolution of the magnetic field structure, the polarization maps can be 
calculated at different time-steps and at any orientation of the modeled galaxy. For the first time a 
self-consistent cosmic-ray electron distribution is used to integrate synchrotron emissivity along the 
line of sight. Finally, our maps are convolved with the given radiotelescope beam. We show that 
it is possible to reconstruct the extended magnetic halo structures of the edge-on galaxies (so called 
X-shaped structures). 

Subject headings: ISM: galactic dynamo — magnetic field — radio polarized emission 



1. INTRODUCTION 

Recent deep observations of polarized radio-continuum 
emission at c entimeter wavelengths of edge-on spi- 
ral galaxies (iTiillmann et al. I l2000t ISoida I 120051 
iHeesen et al. I l2005( ) revealed characteristic X-shaped 
magnetic field structure. Many galaxies showed magnetic 
field oriented in the halo at a large angle to the galaxy 
disk. Earlier observations, mostly limited by sensitiv- 
ity to narrow area along the major axis, demonstrated 
main ly plane-parallel component (e.g. iDumke et al. I 
I1995D . The polarization vectors of NGC 5775 close to 
the galactic plane are parallel to the disk, but in the 
halo they extend up to 10 kpc above the disk form- 
ing so called X-shaped structure in all four quadrants. 
Such picture most probably reflects the qu adrupole con- 
figura tion of magnetic field in NGC 5775 ISoida et al. I 
(|2008f ). The observed degree of polarization is very high 
(locally reaching even up to 50%). It resembles high reg- 
ularity of magnetic field there. Such structures seem to 
be common among edge-on spiral galaxies (NGC 5775, 
NGC 46 66, NGC 4217 NGC 3628, NGC 253, NGC 891 
(see e.g. ' Dahlem et aTl |1997: Sukumar & Allen ' Il991l : 
ISoida 2001). In addition, Tiillmann ct al. (2003) de- 
tected that NGC 5775 shows differential rotation in the 
direction perpendicular to the disk. So far, there is no 
physical explanation for those extended structures. 

An attempt to explain the structure of polarized vec- 
tors by the dipolar dyna mo wave in the NGC 5775 was 
made bv lSokoloff I (|2002t l. He used a solution of the dipo- 
lar classical turbulent dynamo equation and as a result 
he got structures not confirmed by the observations so 
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far ISoida et aTl ((200l . 

The principle of the action of the CR - driven dynamo 
(iParker I Il992t iHanasz fc Lesch I 120001 : iHanasz et aTl 
1200 4?) is based on the cosmic ray (CR) energy supplied 
continuously by supernova (SN) remnants. Due to the 
anisotropic diffusion of cosmic rays along the horizon- 
tal magnetic field lines, cosmic rays tend to accumulate 
within the disc volume. However, the configuration strat- 
ified by the vertical gravity is unstable against the Parker 
instability. Buoyancy effects induce the vertical and hor- 
izontal motions of the fluid and the formation of undu- 
latory patterns - magnetic loops in the frozen-in, pre- 
dominantly horizontal magnetic fields. The presence of 
rotation in galactic disks implies a coherent twisting of 
the loops by means of the Coriolis force leading to the 
generation of the small-scale radial magnetic field compo- 
nents. The next phase is merging the small-scale loops 
by the magnetic reconnection process which forms the 
large scale radial magnetic fields. Finally, the differential 
rotation stretches the radial magnetic field amplifying 
the large-scale azimuthal magnetic field. The coupling 
of amplification processes of the radial and azimuthal 
magnetic field components results in exponential growth 
of the large scale magnetic field with timescales of galac- 
tic rot ation period (140 Myr) as shown by Hanasz ct aL^ 

(pool . 

In this paper we present the results of galactic disk 
modeling which is constructed from a set of local volumes 
placed at different distances from the galactic center. We 
show the results of integration of synchrotron emissivi- 
ties (Stokes parameters I, Q and U) along the line of 
sight (LOS), and compare them to the observed radio 
images of the galaxy NGC 5775. Our results are com- 
pared t o polarized radio contin uum observations made 
at 6 cm (jTiillmann et al. Il2000( ). 

2. MODELS OF THE COSMIC-RAY DRIVEN DYNAMO 
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The first complete 3D numerical mod el of the CR- 
driven dynamo has been demonstrated bv [ Hanasz et al. I 
(|2004 120051 120061) : lOtmianowska-Mazur et al. I (|2007D . 
Our model of the cosmic ray driven dynamo includes 
the following physical elements: the ionized gas and 
magnetic field described by resistive MHD equations, 
the cosmic ray component described by the diffusion- 
advection transport equation (see iHanasz fc Lesch I 
120031 for the details of the numerical algorithm), 
cosmic rays diffusing anisotropi c ally al ong magnetic 
field lines ()Giacalone fc Jokipii I 119991 : iJoki pii et al.^ 
mil), supernova remnants exploding randomly in the 
disk volume, th e finite (curren tly uniform) resistivity 
of th e ISM (see IHanasz et al. I 2002: Hana sz fc Lesch I 
20031: iKowal Hanasz fc Otmianowska-Mazurl 120031 : 
Tanuma et al. I I2003fl and the realistic vertical disk 



gravity and rotation ()Ferriere I I1998D . In the present 
models gas motions arise only from the cosmic-ray 
pressure gradients. We neglect any shock effects, like 
gas heating. It is caused by limitations of currently used 
explicit algorithm for cosmic ray difusion, although we 
plan to incorporate shock effects in our future papers. 
A complementary work incorporating affects of dynamo 
powered directly by supernova-driven turbulence, but 
without taking into accou nt the cosmic rays , has been 
recently published by Gres sel et al. I ()2008al lbl). In the 
present paper we attempt to examine observational 
properties of our model of cosmic ray-driven dynamo. 

The system of coordinates x, y, z corresponds locally to 
the global galactic cylindrical system r,(j), z. The bound- 
ary conditions are periodi c in the Y-direction, sheared in 
the X-direction (foUowing lHawlev et al. |[l995f ) and open 
in the Z-direction for fluid quantities and magnetic field 
components. The boundary conditions for cosmic rays 
are fixed (ecr = 0) on the Z-boundaries. 

We present a new set of numerical models of the CR- 
driven dynamo computed in cubes situated at different 
distances from the center of the galaxy (see Fig. 1). We 
intend to obtain the picture of the whole galaxy from 
these local cubes. In the inner part of the galaxy a vari- 
ety of different physical mechanisms occur (such as cen- 
tral activity, bulge influence etc.), most of which are not 
taken into account by our model. For this reason, we 
start our calculations at 2 kpc radius. The physical pa- 
rameters used in_the calculations are assumed according 
to the paper of Ferricre ( 1998) which refers to the Milky 
Way (see Fig. 1). Our numerical simu lations are per- 
formed using the Zeus-3D MHD code (jStone fc Normani 
Il992airbl ) with the cosmic rays extension made by Hanasz 
fc Lesch (2003). 

The full set of equations describing the model includes 
the set of resistive MHD equatio ns completed by the cos 
mic ray transport equation (see IHanasz et 
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p=(7-l)e, 7 = 5/3 (5) 

where q = — d InfJ/d InR is the shearing parameter, (R 
is the distance to galactic center), rj is the resistivity, 7 
is the adiabatic index of thermal gas, the gradient of cos- 
mic ray pressure Vper is include d in the equation of mo- 
tion (see lBerezinskii et al. Ill990l e.g.) and other symbols 
have their usual meaning. The uniform resis tivity is in- 
clude d only in the induction equation (sec Han asz et al. I 
l2002f) . The thermal, ionized gas component is treated as 
an adiabatic medium. 

The transport of the cosmic ray component is described 
by the diffusion-advection equation 



V {e„V) - V (kVe^r) -Pcr (V • V) + Qsn, (6) 



deer 
dt 

where Qsn represents the source term for the cosmic ray 
energy density: the rate of production of cosmic rays 
injected locally in SN remnants and 



Pc 



(7cr - l)ecr, 7cr = 14/9. 



(7) 



The adiabatic index of the cosmic ray gas jcr ^^nd the 
formula for diffusion tensor 



{K\\- Ki_)n^nj, n^ = B,/B, (8) 

are ad opted following the argumentation by iRvu et al. I 
((20031) . 

In order to construct the polarization maps of an edge- 
on galaxy characterized by a high halo, our local simu- 
lation is extended to 4 kpc above and below the disk. 
The local cube X-size is 500 pc, while Y-size is 1000 pc. 
The cell sizes are equal to dx=dy=dz=20pc. The reso- 
lution of our local calculations is given by 25 x 50 x 400 
grid points. However, we are primarily interested is in 
the large scale magnetic field amplification and structure. 
Because CR gas propagates diffusively and gas dynam- 
ics does not involve shock waves, the currently assumed 
grid resolution is sufficient. This statement is verified by 
resolution studies - simulations performed with the cell 
size (10 pc)'^ up to (20 pc)"^ lead to convergent results. 

The disk rotation was defined by the values of the 
angular velocity at the different galactocentric radii 
ranging from 0.05 Myr""'^ at Rg = 5 kpc down to 
0.025 Myr~^ at Rq = 10 kpc (see. Fig. 1). The value of 
the shearing parameter is (7 = up to a radius of 3 kpc, 
and linearly increasing to g = 1 at 5 kpc, as it is pre- 
sented in Fig. 1. The cosmic ray diffusion coefficients 
assumed in the simulations are: i^n = 3 x lO^^cm^s^^ 
and K± — 3x 10^^ cm^ s~^ . These values are scaled down 
by an order of r nagnitude with resp ect to expected real- 
istic values (e.g. IJokipii et al. 1(199 9) due to the timestep 
limitation in the currently used explicit algorithm for the 
diffusion equation. The assumed value of the resistivity 
coefficient 77 is 3 x 10^^ cm^ s~^. The values of other input 
parameters are depicted in Fig. [1] 

3. CONSTRUCTION OF SYNTHETIC RADIO MAPS 

As the result of our local calculations we obtained 3D 
rectangular cubes of the magnetic field and the cosmic- 
ray energy density. Since the CR component is described 
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Fig. 1. — Radially dependent parameters for all our local calcu- 
lations. Subsequent panels (counting from the top) display angular 
velocity £7, shearing parameter Qsheary column density cr and 
supernova rate /sAf, vs. galactocentric radius Rg, deduced from 
the model by Ferriere (1998). 

as adiabatic relativistic gas instead of the full momentum 
distribution function, the energy density of cosmic rays 
can be considered as proportional to the number den- 
sity of CR nucleons. The number density of cosmic ray 
electrons ne,cr is typically assumed to be of the order of 
1 % of the number density of cosmic ray nucleons n„_cr- 
Therefore rte.cr can be assumed proportional to Ccr- The 
local cubes at given radius and chosen time-step are repli- 
cated into subsequent cylinders, and cylinders are com- 
bined together into the full galactic disk (at the radii 2 
- lOkpc, and of 8kpc thickness). The synthesized disk 
can be oriented according to the inclination and posi- 
tion angle of any real galaxy (in the case of NGC 5775 
it is 86° and 145°, respectively). After computing the 
magnetic field component perpendicular to the line of 
sight B± and relativistic-electron density Ue^cr from the 
simulation data we calculate the synchrotron emissivity 
at each point. Following the standard formula (see e.g. 
iLongair Ill994[ ) we obtain Stokes parameters I, Q, and U: 
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dl 
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6/ 

pti COS 2x cos A — sin A 
pei sin 2x sin A cos A 




(9) 

where the synchrotron emissivity is 

and A denotes the Faraday rotation angle. Integrating 
the Stokes parameters along the line of sight we obtain 
the map of synchrotron emission of the simulated galaxy. 
We set 7 = 2.8 and p = 75%. 

In addition, assuming that the distribution of thermal 
electrons Ue^th is proportional to the gas density, we can 



account for the Faraday rotation of polarized emission 
along the line of sight. The Faraday rotation in small 
distance dl is 



A CX Tie 



,thB\\dl. 



(11) 



Finally, we convolve resulting I, Q, and U maps with 
Gaussian beam of HPBW of 20" and calculate polarized 
intensity and vectors of magnetic polarization {B) for 
direct comparison with real observations of NGC 5775. 

4. RESULTS 

The evolution of the magnetic field energy, the cosmic 
ray energy and the mean values of the magnetic pitch 
angles calculated in the subsequent rings for five chosen 
distances from the galactic center is presented in Fig. [21 

In all local models with nonvanishing differential rota- 
tion ( Qshear 7^ 0) i.e. at the radii i? > 3 kpc, a fast growth 
of the total magnetic field energy is observed (see Fig. [21 

- top panel). Within the ring R = 4-5 kpc the mean 
growth time of the large-scale magnetic field, measured 
in the phase of exponential amplification of the magnetic 
field (since t = 0.2 Gyr till t = 0.9 Gyr) is the shortest 

- about 0.15 Gyr. The growth time gradually increases 
with the radius up to about 0.3 Gyr for the outermost dis- 
tance from the galactic center R = 10 kpc. The fastest 
growth of the magnetic energy at 4.5 kpc coincides with 
the largest SN rate and high shear in comparison with 
other disk areas. For the smallest radii {R < 3 kpc), 
with assumed rigid rotation, the magnetic field is not 
amplified. This justifies our decision of not taking into 
consideration rings close to the center of the galaxy. 

After the exponential amplification phase the magnetic 
field energy enters the saturation phase - the earlier, the 
quicker growth rate was observed. The fastest growth of 
the CR energy (see Fig[51 - middle panel) is present at 
the radius 4.5 kpc (solid line) due to the presence of the 
highest SN rates there, similarly to the magnetic field 
energy behavior at this galactocentric distance. In our 
simulations we do not obtain the energy equipartition 
between the cosmic rays and the magnetic field. The CR 
energy exceeds the magnetic energy by 1 2 orders of 
magnitude. Possible reasons for this excess are: the low 
values of the CR diffusion coefficients, as w ell as, periodic 
boundary conditions (see iHanasz et al~l [2008. for more 
detailed discussion). The bottom graph in Fig. [2l shows 
the time evolution of the mean magnetic pitch angles at 
the same radii as for the magnetic field and the cosmic 
ray energies. In the beginning stage of evolution (until 
0.2 Gyr) all curves exhibit maxima up to 10°. Then, at 
two distances from the galactic center, 4.5 kpc and 6 kpc, 
the values of the mean magnetic pitch angles change sign 
to the negative ones. In the end, after 800 Myr, the pitch 
angles change again sign to positive values and stabilize 
for all curves above 5°. 

Fig. [3l shows the radial dependency of the magnetic 
field energy (top panel), the cosmic rays energy (middle 
panel), and the mean pitch angle values (bottom panel) 
for four chosen time steps. We notice that the magnetic 
field energy grows with time and decreases with the ra- 
dius, as we present in the Figure above. In addition, all 
curves (Fig [3l the middle panel) showing the radial de- 
pendence of the cosmic ray energy exhibit the maxima 
in their central parts (3-4 kpc). The highest increase is 
observed for times between 1.5 and 2 Gyr and for the 
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Fig. 2. — Evolution of the magnetic field energy, the CR en- 
ergy and the mean pitch angle values at five chosen distances (R) 
from the galactic center: 3kpc (the dotted lines), 4.5 kpc (the solid 
lines), 6 kpc (the dashed lines), 8 kpc (the dot-dashed lines), and 
10 kpc (the double dot-dashed lines) 

radius about 4.5 kpc. As we explained above, the CR 
energy maxima are connected with the largest rate of 
SN between 4 and 5 kpc (see Fig.[TJ 

Due to small shear in the central part (see Fig.[T|), the 
mean magnetic pitch angle (Fig. [3l bottom) reaches the 
highest values for all chosen time steps. The smallest or 
negative pitch angle values are obtained at time 0.5 Gyr 
at the radius 4.5 kpc (like in Fig. [2] ~ bottom panel). 
We can conclude that when the dynamo is still under 
development the pitch angles are low or even negative. 
Later on, when the energy grows the pitch angles attains 
about 5° for all radii. 

The main goal of this paper is to compare the radio 
continuum emission obtained from the cosmic ray dy- 
namo model with the observed map of a real galaxy. In 
Fig. |4] we present the face-on view of our model map at 
the final time of 2 Gyr. The figure shows polarization 
vectors superimposed onto the contours and grayplot of 
the polarized intensity. We see that the central part ex- 
hibits the highest pitch angle values. Further out smaller 
pitch angles are observed, in agreement with the bottom 
pannels of Figs. [2] and [3l This fact can be explained byh 
the lower shear in the galactic center than outward in 
the disk. 

Fig. [5] shows calculated maps of the vectors of the po- 
larized emission superimposed onto the isolines and gray- 
plots of the same quantity with orientation on the sky- 
plane as the spiral galaxy NGC 5775 (inclination of 86° 
and position angle of 145°) at selected time steps. One 
can notice that the extended structures of the polariza- 
tion vectors in the modeled galaxy are present from the 
early stages (500 Myr) of the evolution. They appear at 
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Fig. 3. — Radial dependency of the magnetic field energy, the 
cosmic ray energy and the mean pitch angle values for four cho- 
sen time steps: 0.5 Gyr (dotted lines), 1 Gyr (solid lines), 1.5 Gyr 
(dashed lines), and 2 Gyr (dot-dashed lines). 
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Fig. 4. — Face-on map of the polarized emission of the modeled 
galaxy. Contours and gray scale show polarized emission intensity, 
vectors are of directions of apparent magnetic polarization vector, 
and length proportional to the intensity 
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Fig. 5, — Polarization maps of the modeled galaxy at selected 
time steps: 500 Myr (top panel), 1000 Myr (middle panel), and 
1500 Myr (bottom panel) oriented as the real galaxy NGC 5775. 
Contours and gray scale show polarized emission intensity, vec- 
tors are of directions of apparent magnetic polarization vector, and 
length proportional to the intensity 



Fig, 6. — Map of the modeled galaxy at the latest (2Gyr) time- 
step of the evolution (upper panel) oriented as the real galaxy 
NGC 5775 (lower panel). Contours and gray scale show polarized 
emission intensity, vectors are of directions of apparent magnetic 
polarization vector, and length proportional to the intensity 



large distances from the disk forming an X-shaped struc- 
ture, similarly to the observed map of the polarization 
emission of the edge-on galaxies. The most extended 
structures are apparent at the late time steps of our sim- 
ulations 1500 Myr (the bottom panel) and 2000 Myr (see 
Fig. [SI the top panel). 

In Fig. [S] we compare maps of polarized emission for 
NGC 5775 (the lower panel) and synthetic maps for our 
model of this galaxy (the upper panel). In both pictures, 
isolines of polarized intensity are superimposed onto the 
grayplot of the same quantity, together with the vectors 
of directions of the apparent magnetic polarization with 
the length proportional to the intensity. The polariza- 
tion vectors in the central part of the bodies are parallel 
to the disk in both modeled and real galaxies. The ex- 
tended structures, polarized perpendicularly to the disk 
plane, are visible in both maps as well. However, in the 
real NGC 5775 the extensions are separated less than in 
our model. The depolarized canals, clearly seen in the 
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Fig. 7. — Rotation measure plot along the major axis of the 
galaxy — calculated from the model at four time steps and the real 
observations of NGC 5775. Negative distances denotes left (south- 
eastern) side of the galaxy 

model map are barely visible in eastern and western ex- 
tensions in the real galaxy. This suggests that in the real 
galaxy the magnetic field changes its orientation more 
smoothly than in the model. Clear differences between 
the model and real observations appear close to the ma- 
jor axis of the galaxy. Narrow depolarized channel is seen 
on the lower panel of Fig. [6l It can be explained by rela- 
tively strong Faraday effects caused by non-axisymmetric 
magnetic field configuration present in the front side (th e 
south-western side) close to the disk plane ()Soida II2005D . 
Another reason for this behavior could be connected to a 
local enhancement of magnetic field and / or thermal elec- 
tron density due to the infiuence of a spiral arm. Our 
model does not include any non-axisymmetric features. 
Furthermore we do not include the very central part in 
our modeled galaxy. This results in dumping magnetic 
field at the radii i? < 3 kpc. In real galaxies any non-axial 
(vertical) motions can strongly affect magnetic field con- 
figuration. 

We compare our model with observations at relatively 
high frequency (4.86 GHz), where Faraday effects usually 
are small. Nevertheless we include such effects in our cal- 
culation of Stokes parameters (Eq. 9). Calculating model 
maps of polarized emission (intensity and directions) at 
two frequencies we were able to construct the synthetic 
maps of rotation measures (RM) as well. In the RM 
calculations we apply the density of gas taken directly 
from our simulations. A cut of the rotation measure dis- 
tributions along the major axis of both modeled at four 
different time steps and observed galaxy (NGC 5775) for 
comparison are presented in Fig. [71 Our models repro- 
duce the general trend of the rotation measure. The 
deviations - most prominent jumps around 6 kpc of both 
sides coincide with the narrow depolarized canals on the 
polarized intensity map (Fig. [6] - upper panel), where 
the polarization vectors change their direction rapidly 
(that leads to large uncertainties of RM determination) . 
Other differences between the model and observations 
(in the center and at the peripheries) are possibly caused 
by non-axisymmetries in magnetic field and thermal gas 
distribution (such as spiral arms), certainly present in 
NGC 5775 and not included in our model. 

5. DISCUSSION 



The extended structures appear in the modelled galaxy 
as a result of variety of physical processes. First of all, the 
buoyancy driven by CRs transports the magnetic field 
together with CR gas to the halo. We find that the winds 
in halo build the vertical magnetic component. Those 
winds reach speeds up to l OOkm/s) and form q uickly at 
the beginning stage of (see lHanasz et al.1l2004D . 

Our local-box simulations show excessive energy den- 
sity of CRs with respect to to the magnetic energy. Pos- 
sible reasons of this deviation from equipartition are the 
reduced CR parallel and perpendicular diffusion coef- 
ficients and the assumed horizontal periodicity of the 
computational domain. The latter assumption leads to 
trapping of CRs by the horizontal magnetic field, (see 
iHanasz et a"ni2008[ ). 

Nevertheless, for the first time we apply the modeled 
CR distribution to perform calculations of the polarized 
emission in our modeled galaxy. The obtained results 
are in large extent in agreement with observations. The 
model presented in this paper demonstrates that, even 
if the modeled disk starts from 2 kpc from the center of 
the galaxy, the magnetic field still can take the form of 
X-shaped structures in the halo. 

The maps of face-on model show that the pitch angles 
of the polarization vectors are about 5°. This value is a 
bit smaller than the angles usually observed in the spiral 
galaxies. We underline the fact, however, that our model 
does not take into account effects of spiral arms in the 
galaxy. The arms understood as manifestation of density 
waves are known to influence the magnetic field structure 
in galactic disks significantly (Beck & Hoernes 1996[). 

Our model calculates polarized intensity only. It does 
not include any turbulent structures below relatively 
large grid size (20, pc). This results in lack of total 
power emission map comaprable with real observations. 
It does not allow any reasonable discussion of depolar- 
ization other than caused by beam smoothing and/or 
Faraday effects. 

The Faraday effects computed in our model does not 
influence the overall results. Differences between the ro- 
tation measure model and observations are certainly due 
to limitations of axisymmetric model imposed by extend- 
ing the local-box into the global galactic model. 

One could make the final picture more realistic via the 
incorporation of supernova remnants. The incorporation 
of heat output from supernovae may lower the cosmic 
ray excess. We suggest that global simulations of cosmic 
ray-driven dynamo may solve the problem as well, how- 
ever both mentionned solutions of the CR excess problem 
require much larger computational resources and other 
numerical algorithms than those currently available for 
our project. It seems however that the large-scale qual- 
itative features, namely the X-type structures can be 
succesfully reproduced with the present setup. There- 
fore we would prefer to demonstrate the observational 
properties of our present incomplete model of cosmic-ray 
driven dynamo, rather than to mix quantitative solu- 
tions with parametrized effects of supernova remnants. 
In the present simulations with the shearing box approx- 
imation, it is not yet possible to model the vertical shear 
presumably taking plac e in the case of NGC 5775 galaxy 
(jTiillmann et al. il2000[ ). Therefore the global CR-MHD 
simulations of galactic disks free of limitations of the 
shearing box should be addressed in the next step of 
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development of our models. 

Five galaxies NGC 891, NGC 3628, NGC 4217, 
NGC 4 666, and NGC 5775 out of six shown in 
ISoida I ()2005f ) characterized by the X-shaped struc- 
tures of magnetic fi eld c o nfiguration, have flat ro- 
tation curve (iSofu el 1 19961 : iRhee fc van Albadal 119661 : 
iMathewson et al. i ri992: H eald et al. II2007D . One case 
(NGC 4631) exhibits different magnetic field orienta- 
tion with field lines crossing the disk plane. This 
galaxy also is different when the rotation curve is taken 
into consideration, ft rotates rigidly to the large radii 
(jde Vaucouleurs fc de Vaucouleurs 1119631 ). The lack of 
differential rotation favors development of vertical struc- 
tures at all galactic altitudes and reduces the speed of 
transformation of the radial component of magnetic field 
into the azimuthal one Siejkowski (2008) . We suggest 
that the peculiar structure of the polarized emission 
which is observed in NGC 4631, where vectors cross the 
disk plane, can result from the cosmic-ray driven dy- 
namo. This state ment is supp orted by a number of ob- 
servations ^oida_e_£in ()2008D . The galaxy NGC 4631 
drew our attention and we plan to consider it in our fu- 
ture research. 

6. CONCLUSIONS 

In the present paper we compare the maps of polar- 
ized synchrotron radio emission of galaxy NGC 5775 with 
those constructed from the numerical model of the cos- 
mic ray driven dynamo. The main conclusions can be 
summarized as follow: 

• All models of the local cosmic-ray driven dynamo 
computed at the galactic radii between 4kpc and 



lOkpc indicate the fast growth of the magnetic flux 
and the total magnetic energy. 

• The synthetic radio maps of polarized emission 
computed on the basis of our local models exhibit 
vertical magnetic field structures similar to those 
observed in numerous edge-on galaxies. 

• The polarization vectors in the disk plane (face-on) 
form a spiral pattern with the pitch angles about 
5°. The pitch angles from the model are slightly 
smaller than we normally observe in galaxies. How- 
ever, the pitch angle of the mean magnetic field de- 
pends to some extent on the actual parameters like 
e.g. the magnitude of the CR diffusion coefficients. 

• As it was expected for observation at relatively high 
frequencies, the Faraday effects considered in our 
model have no significant importance on the final 
results. 

• The major simplification of our model is caused 
by assumed axial symmetry, which is necessary for 
construction of the global model from local-box 
simulations. 
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